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Structural instabilities and sequence of phase transitions in SrBi 2 Nb 209 and 
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Despite their structural similarities, SrBi2Ta209 (SET) and SrBi2Nb209 (SBN) undergo a dif¬ 
ferent sequence of phase transitions. The phase diagram of SBT as a function of the temperature 
includes an intermediate phase between the high-temperature phase and the ferroelectric ground 
state, while in the Niobium compound the intermediate phase is suppressed and a single transition 
between the high- and low temperature structures is observed. We present ab initio calculations that 
reveal the relevance of a trilinear coupling between three symmetry-adapted modes to stabilize the 
ground sate in both compounds, being this coupling much stronger in SBN. Within the framework 
of the phenomenological Landau theory, it is shown that by solely increasing the strength of the 
trilinear coupling the topology of the phase diagram of SBT can change up to suppress the inter¬ 
mediate phase. Monte Carlo simulations on an idealized Hamiltonian confirm that the trilinear 
coupling is the key parameter that determines the sequence of phase transitions and that for higher 
dimensionality of the order parameters the stability region of the intermediate phase is narrower. 

PACS numbers: 77.80.-e, 77.84.-s, 63.20.Ry, 71.15.Mb 


I. INTRODUCTION 

The Aurivillius phases are layered bismuth compounds 
that obey the general formula Bi2mA„_mB„03„^^— . The 
family includes many members that present ferroelectric- 
ity at room temperature, and they have been widely 
studied for potential technological applications mainly 
in thin film nonvolatile memoriesi^^i SrBi2Ta209 (SBT) 
and SrBi2Nb2 09 (SBN) are members of the family, with 
n = 2 and m = 1 and their structure is formed by al¬ 
ternating two SrMOa (M=Ta, Nb) perovskite blocks and 
one Bi202 slab (Fig. [T|). 

Tantalum and Niobium present very similar physi¬ 
cal and chemical properties, including the valence and 
atomic radii. As a result, SBT and SBN present 
isomorphous polar structures at room temperature^^iS 
and analogous mechanical and electrical characteristics. 
However, their phase diagrams are qualitatively different: 
on increasing the temperature SBT undergoes a phase 
transition to a non-polar orthorhombic phase that does 
not arise in SBNi^^ii 

Previous first-principles calculations by Perez-Mato et 
al revealed that the trilinear coupling of two primary 
unstable modes and a secondary hard mode is critical to 
stabilize the ground state of SBT.™ This work provided 
a simple and plausible scheme in contrast with others 
alternatives where the eventual existance of negative bi¬ 
quadratic coupling triggers the simultaneous condensa¬ 
tion of several modesJ^ Later, the key role of the trilinear 
coupling has been found in several other compoundsi^^SS 
and its influence has been proposed as the origen of the 
so-called avalanched phase transitions and the hybrid 
improper ferroelectricityi^ 

In this work we revisit the SBT and extend the ab ini¬ 
tio study to SBN in order to investigate the role of the 
trilinear coupling in the stabilization of the ferroelectric 
structure and its influence in the sequences of phase dia- 



FIG. 1: Structure of SBT and SBN in the tetragonal phase 
(M=Ta,Nb). 


grams. 


II. SYMMETRY CONSIDERATIONS, 
DISTORTIONS AND PHASE TRANSITIONS 

The high- and low-temperatures structures of SBT 
and SBN are isomorphous: tetragonal with space group 
lA/mmm (No. 139) at high enough temperatures^ 
and orthorhombic A2iam (No. 36) at low tempera¬ 
tures. In the case of SBN there is a direct phase tran¬ 
sition from the tetragonal to the ferroelectric structure 
at 800K;2ii£ whereas SBT presents an intermediate non- 
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polar orthorhombic Amam phase (No. 63) between 500- 
600K and TTOKA The cell parameters of the two or¬ 
thorhombic phases are related to the tetragonal ones ac¬ 
cording to Co = ct, ao = ax -b bx and bo = ax — bx, 
which means that the volume of the unit cell is doubled 
in the tetragonal-to-orthorhombic transitions. 

The tetragonal space group {lA/mmm) is a supergroup 
of the two orthorhombic space groups, the intermediate 
phase of SBT (Amam) and the ground state of both 
compounds {A2iam). In consequence, the low tempera¬ 
ture structures can be described in terms of symmetry- 
adapted distortions with respect to the parent tetragonal 
structure;^ TableUshows this decomposition in terms of 
the symmetry-adapted modes that break the tetragonal 
symmetry for the experimental structures of SBT and 
SBN obtained with the aid of the tool AMPLIMODES 
(www.cryst.ehu.es).— The hierarchy of the modes is simi¬ 
lar in both compounds: the amplitude of is dominant 
and can be considered as the primary order parameter 
and the distortion seems to play a secondary role. 
The analysis of the T^ polar distortion is more delicate 
because due to the polar character of the ground state 
its amplitude depends on the origin. In Table |T] the cho¬ 
sen origin does not change the arithmetic center of the 
parent structure. According to the isotropy subgroups 
the intermediate phase of SBT is naturally explained as 
a condensation of the mode, and in terms of the Lan¬ 
dau theory of phase transitions the lA/mmm —^ Amam 
transition can be continuous. The second transition to 
the ferroelectric structure should then consist on the con¬ 
densation of the polar mode. In the case of SBN a direct 
transition between the tetragonal and the ferroelectric 
phase can be solely described by a simultaneous conden¬ 
sation of at least any pair of the three modes, and accord¬ 
ing to Landau theory, it should be first order: the three 
relevant modes condense simultaneously in an avalanche 
phase transitional. 

The orthorhombic strain is usually very small in the 
Aurivillius compounds and SBT is not an exception. 
The experimental values of cell parameters in SBT are 
ar = br = 3.917 A and cx = 25.114 A for the tetrag¬ 
onal structure^ and ao = 5.531 k , bo = 5.534 A and 
Co = 24.984 A for the ferroelectric phasei^ The deforma¬ 
tion of the cell is very small with ^plaTjcT ^ aojco ^ 
bo I Co (all within a 0.5 % margin), and the influence of 
strain can be considered as a residual effect. Up to our 
knowledge, there is not such an experimental data avail¬ 
able for the tetragonal cell of SBN, and the orthorhombic 
cell parameters are ao = 5.519 A and bo = 5.515 A and 
CO = 25.112 A4 

A complete description of the atomic displacements 
associated with the three modes can be found in Ref. 
fl^ . In summary, among the 18 degrees of freedom that 
constitute a generic distortion that links the tetragonal 
and ferroelectric structure, seven correspond to X ^, and 
they mainly account for tiltings of the oxygen octahe- 
dra around the (1, 1,0 )t direction of the tetragonal cell 
and displacements of the Bi cations along the (1, —1, 0 )t 


Irrep Direction 

Isotropy 

subgroup 


Amplitude (A) 



SBT 


SBN 

exp.— 

6 

exp.— 

calc. 

exp. 

calc. 

Xg" (a, -a) 

Amam{63) 

0.90 

0.86 

1.12 

0.81 

1.14 



(0.99) 

(0.99) 


(0.99) 


Tg (a, a) 

F2mm{46) 

0.51 

0.63 

1.01 

0.74 

1.03 



(0.94) 

(0.98) 


(0.81) 


X+ (a, -a) 

Abam{6A) 

0.26 

0.39 

0.56 

0.63 

0.75 



(0.84) 

(0.99) 


(0.99) 



TABLE I: Summary of the mode decomposition of the exper¬ 
imental and calculated distortions of the ferroelectric phases 
of SBT and SBN with respect to their corresponding par¬ 
ent (calculated and experimental) structures. Amplitudes are 
given for modes normalized within the primitive unit cell of 
the polar structure. The second row shows in parentheses the 
scalar product of the calculated and experimental polarization 
vectors after normalization. 


direction. The octahedra behave as rigid units and, in 
consequence, neighboring octahedra show antiphase tilt¬ 
ings. The polar T^ distortion involves essentially an an¬ 
tiphase displacement of the Bi atoms and the perovskite 
blocks, along the (1,1, 0 )t direction. Among the eight 
independent displacements associated with the dis¬ 
tortion the main amplitudes correspond to tiltings of the 
oxygen octahedra around the z axis and antiphase dis¬ 
placements of the oxygens in the Bi 202 slabs along the 
(1, 1 , 0)t direction. 


III. COMPUTATIONAL DETAILS 

The WIEN2k code^^, based on the full potential 
LAPW-l-lo method, was employed for the ab initio cal¬ 
culations. Exchange and correlation effects were treated 
within the GGA approximation with the Perdew-Burke- 
Ernzehof parameterization^^. The RK^ax parameter, 
which is related to the number of radial basis functions 
used to describe the the inside spheres, was chosen to 
be 7.5 for both compounds. Calculations in the tetrag¬ 
onal symmetry were performed using a Monkhorst-Pack 
k point mesh of 8 x 8 x 8, which is equivalent to 56 inde¬ 
pendent k points in the irreducible Brillouin wedge. For 
calculations in the orthorhombic symmetry a fc-mesh of 
5x5x5, representing 27 independent k points in the 
IBZ, was chosen in order to keep the k point density as 
constant as possible. In the case of SBT the radii of the 
atomic spheres chosen for the calculations were 2.0(Sr), 
2.26(Bi), 1.8(Ta), 1.6(0) bohrs in the tetragonal struc¬ 
ture and 2.25(Sr), 2.3(Bi), 1.8(Ta), 1.6(0) bohrs in the 
orthorhombic A2iam phase. For SBN, the radii were 
2.26(Sr), 2.26(Bi), 1.8(Nb), 1.66(0) bohrs for both the 
tetragonal and the orthorhombic basis. The choice of 
the parameters was preceded by energy difference con¬ 
vergence tests which confirmed their validity. The con- 
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Compound 

Abs. minimum 

^ 3 - 

Ts" 

A 2 + 

SBT 

15.62 

11.72 

5.24 

1.52 

SBN 

18.34 

10.40 

3.68 

3.52 


TABLE II: Depths of the energy wells per formula unit (in 
mRy) of the ferroelectric ground state (first column) and the 
pure distortions and with respect to the relax;ed 

tetragonal parent conhguration. 


vergence criterion for the SCF calculations was of 0.0001 
Ry for energy and 0.1 mRy/bohr for forces. 

First, relaxations of the parent structures and ground 
states for the two compounds were carried out. In 
the case of SET forces remained below 0.02 mRy/bohr 
and 0.15mRy/bohr in the tetragonal and orthorhombic 
phases respectively. For SEN, the tetragonal relaxation 
accepted an accuracy under 0.1 mRy/bohr and all forces 
present in the orthorhombic relaxed structure were under 
0.17mRy/bohr. 

As mentioned in the previous section the orthorhom¬ 
bic strain is very small in the Aurivillius compounds. In 
consequence its effect was neglected and all the calcula¬ 
tions were performed assuming a fixed cell with tetrag¬ 
onal metric. The c parameter of the idealized cell was 
fixed to the experimental value of the ferroelectric phase 
and the a and b lattice constants were forced to be equal 
and preserve the experimental volume of the same phase. 


IV. ENERGETICS AND COUPLINGS 

First, the atomic positions of both compounds were 
relaxed under tetragonal and orthorhombic A2iam sym¬ 
metries, and the distortions of the low-temperature struc¬ 
ture were analyzed in terms of symmetry modes. As 
shown in Table HI the calculated amplitudes are system¬ 
atically larger than the experimental values, what could 
be related to some degree of disorder in the experimental 
structures that is not present in the calculations. Ta¬ 
ble H] also shows the the scalar product of the calcu¬ 
lated and experimental polarization vectors after normal¬ 
ization. The extremely good agreement indicates that, 
apart from a global amplitude, the calculated and exper¬ 
imental relative atomic displacements are essentially the 
same. 

The energy gains of the polar phase with respect to 
the tetragonal one, together with the depth of the en¬ 
ergy wells associated with the pure distortions are listed 
in Table HVl The primary character of the order pa¬ 
rameter is manifest in both compounds. However, the 
role of the other two distortions seems to be different in 
both compounds. In SET, F^ presents a well defined en¬ 
ergy minimun, while the mode is slightly soft at low 
amplitudes and rapidly hardens for medium amplitudes 
(Fig. [2|) reinforcing its secondary character. In SEN, the 
energy wells for the F^ and modes are very similar 


and a clear hierarchy between both distortions cannot be 
established (Fig. [5]). 

The energy variations around the tetragonal configu¬ 
ration in terms of the amplitudes of the three relevant 
modes can be described by a polynomial composed of 
symmetry invariant terms. The knowledge of the coeffi¬ 
cients in the energy expansion allows the quantification 
of the energy for a general distortion and the strength 
of the couplings between modes. The energetic contribu¬ 
tion of a general combination of A^ , F/" and A/^ can 
be expressed up to fourth order by: 

AE = E^-+ Ej.- + E^++ Ey-j.- 

-A3 i g .A2 A.3 i g 

+^x-x+ + x+ + ^x- r- x+ (1) 


The energy due to the pure modes is given by: 


^x. -+ I^xzQ\ 


1 

1 




^x+ - 2 ^x+Qx+ 


The biquadratic couplings read: 

^xrrr ~ ^xrrrQx~Qr- 


^x-x+ - ^x-x+Qx-Qx+ 
^r-x+ = ^r-x+Ql-Qx+ 


and the trilinear coupling is: 


^x-T-x+ - 1 x-t-x+Qx-QtiQx+ 

The coefficients of the polynomials for both compounds 
were determined by fitting the ab initio energies calcu¬ 
lated in more than 60 points of the configuration space 
for each compound (Table Hill) . The most remarkable dif¬ 
ference between both compounds is the magnitude of the 
trilinear coupling, its value in SEN being much larger 
than in SET. The magnitudes of the rest of coefficients 
are quite similar although the stiffness constants of F^ 
and modes approximately exchange their values and, 
in consequence, their role as the least unstable distortion. 

The condensation of two of the three relevant sym¬ 
metry modes is enough to lower the tetragonal symme¬ 
try down to the ferroelectric space group. However, in 
both compounds the simultaneous freezing of any pair of 
modes is penalized by the strong and positive biquadratic 
couplings 5. FiglHKa) shows the renormalization of the en¬ 
ergies due to the biquadratic coupling for SET: the pres¬ 
ence of a nonzero F/" distortion stabilizes the soft 
mode and a simultaneous condensation of both modes 
becomes energetically unfavourable. The inclusion of the 
a/" mode is essential to explain the polar ground state 
of SET. In figEIb) is shown the energy of SET as a func¬ 
tion of the amplitude of the pure Xj mode and when 
a mixed distortion with Qy-=0.90 A and Qp-=0.84A 

^3 ^ 5 
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SBT 

SBN 

(mRy/bohr^) 

-44.26 

-40.42 

Kp- (mRy/bohr^) 

-21.41 

-17.28 

(mRy/bohr^) 

-16.39 

-23.49 

(mRy/bohr^) 

11.03 

10.1 

/3p- (mRy/bohr**) 

5.62 

5.40 

/3y+ (mRy/bohr"^) 

11.50 

10.25 

(5j^-p- (mRy/bohr'*) 

6.20 

7.76 

Sx-x+ (mRy/bohr"^) 

10.60 

12.60 

'^x+r“ (mRy/bohr'*) 

-2.30 

3.83 

(mRy/bohr®) 

^3 ^5 '^2 

-13.81 

-22.55 


TABLE III: Polynomial coefficients of the energy expansion 
of SBT and SBN obtained by least-square fits. 

is applied; the slightly hard mode becomes strongly 
unstable when the two primary modes condense. This is 
a manifestation of the critical relevance of the trilinear 
coupling between the three modes in the stabilization of 
the ground state. The analysis of the couplings in the 
case of SBN yields the same qualitative picture. 

Analogous conclusions were obtained in Ref. [H for the 
case of SBT, i.e., the trilinear coupling between three 
modes is determinant in the stabilization of the ferro¬ 
electric structure. However, the authors studied the cou¬ 
plings between modes that correspond to the eigenvectors 
of the force constant matrix. Particularly, the displace¬ 
ment pattern associated with the static distortion of 
this work involves contributions from two phonons, one 
hard and one soft. In consequence, the values of the cou¬ 
pling constants in Table Hill and some details of Figures 
[2] and [3] cannot be directly compared with the results of 
Ref. I12I . A similar scenario has been found in several 
other compounds,— suggesting that the critical role 
of trilinear couplings to stabilize the ground state could 
be rather common. 


V. APPROXIMATE PHASE DIAGRAM 

The expansion of the energy in the previous section can 
be considered as the zero temperature free energy of the 
system, and a phenomenological free energy can be ap¬ 
proximated under the assumptions of the Landau theory 
of phase transitions. Then, the temperature renormal¬ 
ization is solely contained in the quadratic terms and 
the stiffness constants present a linear dependence on 
temperature, such that Ki = ai{T — Toy) for a given i 
mode, where a and are constants. Thus, the ab ini¬ 
tio stiffness constants are related to the transition tem¬ 
peratures and renormalization constants of the modes by 
Ki = In a first approximation the temperature 

dependence of higher order terms can be neglected and 


the free energy reads: 

F = J2{^iT-To.,)Q^,+P.Qt) 

i 

-f Ev —p— “b Ev — v+ F/p— -j- Ev —p— v+ 

Fig. Sla) and (b) show the phase diagrams of the two 
compounds in the space of the two more negative stiff¬ 
ness constants, k^- and /Cp- for SBT, and k^- and k^+ 
for SBN. The diagrams correspond to a particular sec¬ 
tion where the difference between the stiffness constants 
Kp- and K^+ remains constant, which implies a similar 
temperature renormalization in both order parameters. 
The stability regions of Fig. 2] correspond to the parent 
phase of symmetry lA/mmm (no modes are frozen), the 
A 2 iaTO ground state (the three modes are frozen), and 
three phases associated with the condensation of a single 
mode: Amam for , F2mmm for F^ and Abam for 
X+. 

The phase diagram of SBN presents a finite gate that 
joins the stability regions of the tetragonal and ferroelec¬ 
tric phases. Thus, an appropriate renormalization of the 
order parameters could drive the system directly from 
the ferroelectric state to the parent structure without an 
intermediate phase. However, the topology of the phase 
diagram for SBT forbids this possibility and the inter¬ 
mediate phase must be present in a finite range of tem¬ 
peratures. This result could explain the existence of an 
avalanche phase transition only in the case of SBN. 

The fundamental influence of the strength of the tri¬ 
linear coupling in the features of the phase diagrams 
is shown in Fig EIc). It corresponds to SBT but 
with the trilinear coupling constant changed to that of 
SBN, that is, 7 = —22.55 mRy/bohr^ instead of 7 = 
— 13.81 mRy/bohr^. The higher value of the coupling 
induces a wide direct passage that allows an avalanche 
phase transition. This suggests that the trilinear cou¬ 
pling does not only stabilize the ground states of the two 
compounds, but it also governs their sequence of phase 
transitions. 


VI. DIMENSIONALITY OF THE ORDER 
PARAMETERS AND TRILINEAR COUPLING 

The expressions of the previous section are not com¬ 
pletely general since they do not take into account the 
correct dimensionality of the the order parameters. As 
usual in such calculations, the expansion of the energy in 
terms of distortions from the high-symmetry configura¬ 
tion has been done by fixing the orientation of the multi¬ 
dimensional order parameters along a specific direction, 
the direction that corresponds to one of the observed do¬ 
mains of the low-temperature structure. In the case of 
SBT and SBN the three order parameters, ^ F^ and 
X^, are two-dimensional and the experimental distortion 
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SBT 



SBN 



QW 


FIG. 2: Energy per formula unit relative to the tetragonal structure of SBT (left) and SBN (right) in terms of the amplitudes 
of the and X 2 distortions. 




FIG. 3: Energy per formula unit of SBT along selected lines of the configuration space: (a) in terms of the amplitudes of the 
the pure polar mode (squares), and the same mode after the X^ distortion has been frozen to 1 A (circles), b) in terms of the 
amplitude of the pure X^ mode (squares) and the same mode after the amplitudes of X^ and F^ distortions have been fixed 
to 0.90 A and 0.84 A respectively (circles). The plots have been shifted vertically to fix a common origin. 


corresponds to the specific direction indicated in Table U 
and the symmetry-equivalent ones. 

A general distortions must be expressed as two- 
dimensional vectors and, in order to simplify the nota¬ 
tion, we will use the following correspondences in carte¬ 


sian and polar coordinates: 

Qx- = if = (ipi,(p2) = iPcpCOs0^,p^sm9!^) 
Qp- = rp = {'tjji,^2) = {pijjCOsO^, p^sinO^) 
Qx+ = V = (vi,V 2 ) = {p-n cos Or, , Pn sin Or,) 
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SBT SBN SBT with 7 sbn 




FIG. 4: Phase diagrams for SBT (left) and SBN (center). The variables of the abscises correspond to the second more unstable 
mode, for SBT and for SBN. The diagram for SBN presents a narrow gate that allows a direct transition from the 
ground state to the parent phase. The right panel corresponds to SBT but with the value of the trilinear coupling switched 
to that of SBN and it presents a wide passage that opens the possibility of a direct transition. The arrows show the same 
hypothetical paths along the three phase diagrams. 


The invariant polynomials for general directions of 
the order parameters obtained with the aid of the 
INVARIANTS^! utility are given in Table IIVI in carte¬ 
sian and polar coordinates. The lowest order term that 
includes anisotropy is the third order invariant related 
to the trilinear coupling; moreover, the presence of this 
term is enough to explain the observed ground state. 
For a negative trilinear coupling the absolute minimum 
corresponds to the experimentally observed directions of 
the order parameters (Table |T|, whereas a positive co¬ 
efficient should give a different structure with the same 
space group but Ag", T^ and oriented along (a, —a), 
(a, a) and (—a, a) respectively. 

Following the approach of Ref. [l^, a simple micro¬ 
scopic Hamiltonian that retains the main features of the 
SBT/SBN system can be developed from an extended 
version of the two-dimensional model. The hamilto- 
nian for a single isolated mode can be written as: 

-1]2 + - </.(j)]^ 

* (ij> 

with (f) = (p,ip,ri. Erf, is the depth of the double wells at 
sites i, Crj, are the harmonic couplings between nearest 
neighbours, and the ratio between both determines the 
behaviour of the system from order-disorder {E » C) 
to displacive {E « C). 

The trilinear coupling can be included as an on-site 
interaction: 

= ^^[:pi(z)V’i(f)»7i(*)+ <P2(f)^/’i(j)»72(*) - 

^ i 

and its contribution at the absolute minimum is 
~\l\PvPi’Pv in. the one-dimensional case. In order to 


Cart. coor. 

2 nd order 


3rd order 

Ifllpiril + ‘P 2 lpiri 2 - 'P 2 'ip 2 'ni - V5 iV’2??2 

4th order 

HW \ri\\ 

^>1+ P2, V’l + Vl + Vt, 

Pl<P2'(pl'tp2, ^1^21)1^2 , tpl^2rilP2 

Polar coor. 

2 nd order 

2 2 2 

Ptp! Pt/) ? Pr] 

3rd order 

p<pP4,Prilcos9,p cos{9,f, + 9r,) — s\n9,p sm{9,p — 9r,)] 

4th order 

4 4 4 

P<^i Pi/j ? P?7 ? 

p^(3 -1- cos40^), p^(3 -1- cos40,/,), p^(3 -1- cos40,,), 

2 2 2 2 2 2 

Pv’P'0’ P<pPi7i Pi/iPt?? 

p%p% sin 20,^ sin 20,/,, p%p^ sin 20,^ sin 20,,, 
sin 20,/, sin 20,, 


TABLE IV: Invariant polynomials up to fourth order of the 
components of the order parameters in cartesian and polar 
coordinates, (p, ip and rj correspond to A^, T^ and A^ re¬ 
spectively. 

focus on the effect of the trilinear coupling term and to 
maintain the analogy with the one-dimensional case of 
Ref. [^ . we have not included in the Hamiltonian the 
rest of the fourth order terms that are allowed by sym¬ 
metry. The complete Hamiltonian reads: 

H = Hrp ~\~ Hrf, -f Ffr/ “f Erp^rfi,r] 

The analysis of all the possible combinations of the 
parameters of the Hamiltonian is beyond the scope 
of this work and we have limited the study to some 
particular cases. First, we have considered the same 
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FIG. 5: Amplitudes of the order parameters ip (solid circles) 
and Ip (open circles) for E/C = 10 and three values of the 
trilinear coupling: a) 7 = O.lEy,, b) 7 = 0.2E^ and c) 7 = 
O.SEif,. The amplitudes of rj are the same as those of p) and 
are not shown in the figure. 


displacive/order-disorder degree for the three order pa¬ 
rameters {E^p/C^ = E^/C^ = Erj/Crj = C/E) and 
two cases: order-disorder (E/C = 10) and displacive 
(E/C = 0.1). The dominance of the instability of the 
X/^ (= if) distortion in both compounds has been taken 
into account by assuming = Crj = O.SCip, and in 
consequence E^ = Ej^ = Q.SE^. This choice yields the 
same statistical behaviour for ip and 77 and the relation 
= Trj = O.STip for the transition temperatures of the 
pure modes, i.e., without the trilinear coupling. 

We have studied the statistical properties of the model 
by Monte Carlo simulations. Calculations were carried 
out using cubic supercells containing 24 x 24 x 24 sites 
and statistics were collected during at least 10® Monte 
Carlo steps after a proper equilibration. Fig. [5] shows 
the equilibrium amplitudes of the order parameters for 
the order-disorder case (E/C = 10) and three different 
values of the trilinear coupling. Temperature is given in 
units of Cp. Three different regimes can be observed, 
for low values of 7 [Fig. |S] (a)] both transitions are con¬ 
tinuous and the range of stability of the intermediate 
phase is quite wide; for intermediate values of the trilin¬ 
ear coupling [Fig. [5| (b)] the intermediate phase is stable 
in a narrow range of temperatures and the low temper¬ 
ature transition is first-order, and finally, for stronger 
couplings [Fig. [5] (c)] the intermediate phase disappears 
and a direct discontinuous avalanche phase transition be¬ 
tween the tetragonal and ferroelectric phase occurs with 
the simultaneous condensation of the three distortions. 

Although results are qualitatively similar to those of 
previous works with one-dimensional order parameters, 
the different calculation methods, mean field approxima- 


Compound 

Cas Mn 2 07 

A 2 +( 2 ) 

A 3 -(2) 

r5-(2) 

PbZrOs 

A5-(3) 

E 2 ( 12 ) 

S'2(12) 

Sr 2 MW 06 

rj(3) 

A3+(3) 

A5+(6) 

SBT/SBN 

A 3 -(2) 

r5-(2) 

A 2 +( 2 ) 


TABLE V: Three examples of the multidimensionality of or¬ 
der parameters that couple trilinearly (the labels of the irre¬ 
ducible representations and their dimensions are listed). The 
M cation in Sr 2 MW 06 corresponds to Zn, Ca and Mg 


tion and Monte Carlo simulations, prevent a quantitative 
evaluation of the role of the dimensionality in the stability 
of the phases. Therefore, we have also performed Monte 
Carlo simulations of the one-dimensional model and the 
comparison between the phase diagrams for both cases 
are shown in Fig. [B] An accurate estimation of the tran¬ 
sition temperatures was obtained by inspecting the prob¬ 
ability distributions of the amplitudes of the distortions 
with simulations close to the critical regions. In the dis¬ 
placive regime the intermediate phase is suppressed for 
7 ^ 2.5Ep and 7 ~ 2.0Ep for one and two-dimensional 
order parameters respectively [Fig. ID (a)]. The influence 
of the dimensionality is more remarkable in the order- 
disorder regime where the avalanche phase transition ap¬ 
pears for 7 ^ QAbE^p in the Id case, and 7 ^ 0.3Ep when 
the multidimensionality of the order parameters is con¬ 
sidered. It can be concluded that for higher dimensional¬ 
ities of the irreducible representations associated to the 
order parameters the trilinear coupling is more efficient 
to favour a direct phase transition between the high and 
low-temperature phases. 

Multidimensional order parameters are ubiquitous in 
structural phase transitions. For instance, in completely 
different compounds as the ferroelectric Ca 3 Mn 207 , an- 
tiferroelectric PbZrOa (PZO) and the double perovskites 
Sr 2 MW 06 (M =Zn, Ca and Mg), that do not present 
any dielectric anomaly, the trilinear coupling is critical to 
stabilize the ground state as in SET and SEN. Table IVl 
shows the relevant irreducible representations and their 
dimensions. Whereas the case of Ca 3 Mn 207 is similar 
to the SET/SEN compounds, the high dimensionality of 
the irreps in the case of PZO could favour the observed 
strongly-discontinuous transition observed upon cooling 
from the cubic phase. 


VII. CONCLUSIONS 

Despite the different sequence of phase transitions that 
SET and SEN go through, first principles calculations 
shows that the nature and main features of the instabil¬ 
ities in both compounds are qualitatively very similar, 
and that most of the conclusions about SET obtained in 
Ref. [13 can be extended to SEN. Among the three sym¬ 
metry adapted distortions (A;^, F^ and A^) that take 
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FIG. 6 : Transition temperatures of the one- and two-dimensional order parameters in terms of 7 for a) order-disorder and b) 
displacive systems obtained by Monte Carlo simulations. The grey area shows the stability region of the intermediate phase. 


part to drive the tetragonal lA/mmm parent phase to the 
polar A2iam structure, is dominant with the highest 
amplitude and the deepest energy well (^11 mRy per 
formula unit). The depths of the secondary distortions 
and X 2 are similar in SEN 3.5 mRy p.f.u.), while 
in SET X 2 is hardly unstable and slightly deeper 
(5.2 mRy p.f.u.). 

In both compounds , the simultaneous condensation 
of the primary instability X^ with any of the other two 
distortions is penalized energetically by positive strong 
biquadratic couplings. Thus, the trilinear coupling is the 
key ingredient that allows the simultaneous condensation 
of the three modes and its role is essential to stabilize the 
observed ground state. 

The most noticeable difference between SET and SEN 
is associated with the magnitude of the trilinear term. 
The polynomial expansion of the energy of both com¬ 
pounds shows that the trilinear coupling is much stronger 
in SEN than in SET, suggesting that its magnitude could 
be crucial to suppress the intermediate phase in SEN. 
Moreover, the analysis of phenomenological phase dia¬ 
grams for SET and SEN shows that a higher value of 
this constant is enough to change the topology of the 
phase diagram, allowing a direct phase transition from 
the parent phase to the ground state. 


Finally, a simplified cj)'^ Hamiltonian that retains the 
symmetry requirements and the correct dimensionality 
of the order parameters has been developed for the 
SET/SEN system. Monte Carlo simulations do not 
show qualitatively differences in comparison with the 
one-dimensional case^^: the different sequence of phase 
transitions observed in both compounds can be repro¬ 
duced by changing the strength of the trilinear coupling. 
However, according to the present calculations, the in¬ 
crement of the fluctuations associated with the higher 
dimensionality of the order parameters tends to favour 
the suppression of the intermediate phase and to reinforce 
the first order character of the direct transition between 
the high- and low-temperature structures. 
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